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Abstract 



We investigate the BCS and generalized BCS theories in the relativistic 
quantum field theory. We select the gauge freedom as U(l), and introduce a 
BCS-type effective attractive interaction. After introducing the Gor'kov for- 
malism and performing the group theoretical consideration of the mean fields, 
we solve the relativistic Gor'kov equation and obtain the Green's functions in 
analytical forms. We obtain various types of gap equations. 

I. INTRODUCTION 

The most important turning point in the history of the theory of superconductivity is 
the appearance of the BCS ( Bardeen-Cooper-Schriefer ) theory [1]. It is in fact the first 
success giving an explanation of the superconductivity from the microscopic dynamics of 
the interacting many-fermion system. From a view point of the condensed matter theory, 
the BCS theory can be interpreted as a simplification of the electron-phonon theory, which 
reveals the fundamental origin of the mechanism of the generation of the superconducting 
state. However, the influence of the BCS theory is much more wide. It is the famous fact that 
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its influence reached not only to the condensed matter physics, but also to the elementary 
particle physics and the nuclear physics [2~6]. 

The BCS theory succeeded in explaining the superconductivity in metals and alloys. But 
in the development of the high-T c and heavy-fermion superconductivities, the BCS theory 
could not cover them completely. This situation caused new discussions on the mechanism 
of the origin of these superconductivities [7] . Apart from these discussions on the mechanism 
of the origin, various efforts also have been made to extend ( generalized ) the BCS theory 
to understand these systems [8~15] . In the generalized BCS theory, the question of the 
mechanism of the origin is set aside for a while, and the emphasis is placed on the symmetries 
of the pairing interaction and the mean fields. They extended the treatment of the Cooper 
pairs, not only to the singlet pairing but also to the triplet pairing, in which the mean 
fields have the vectorial nature; they have finite angular momenta. The difference of the 
attitudes between the theories of the mechanisms of superconductivities and the generalized 
BCS approach has an important meaning in the theory of superconductivity. 

Until now, almost all the theories of condensed matter physics are constructed on the 
nonrelativistic quantum theory. On the other hand, in the electronic structure calculations 
based on the band theory, it is obvious that, in the cases of heavy elements, we have to 
take into account several relativistic effects for correct descriptions even qualitatively. For 
example, even in the case of slow electrons like ones with a kinetic energy of 100 eV ( 
the electron mass is about 0.511 MeV ), we need the 4-component Dirac equation for an 
qualitatively correct description of the system [16]. Needless to say, the targets of condensed 
matter theory are more or less to include the relativistic effect. However, the study on the 
relativistic effects in condensed matter is not so much performed, except the considerations 
of the spin-orbit coupling [17]. 

Under this situation, recently papers which assert the necessity of the relativistic treat- 
ment for superconductivity in condensed matter were published by K. Capelle and E. K. U. 
Gross [18~24]. The heavy elements like Au, Pb, Ir or I become superconducting states [25], 
and simultaneously they are under strong relativistic effects. The high-T c or heavy-fermion 
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superconductors normally contain heavy elements like La, Ce and U. Especially in the 
heavy-fermion superconductors, the f-electrons which are under strong relativistic effects, 
become the superconducting state. In principle, the description of the electronic structure 
for these systems requires the relativistic treatment. It is also a well known fact that, the 
spin-orbit effect plays a role in the NMR spectra ( the Knight shift ) and also in the impurity 
effects in superconductors [17]. Naturally it becomes interesting to treat these physics with 
the strong relativistic effects, beyond making a small corrections. There are no reasons that 
those systems which have strong relativistic effects do not show relativistically distinct su- 
perconductivity than the nonrelativistic one. In the works of Capelle and Gross [18,19], the 
necessity of the relativistic treatment of the theory of superconductivity in solid state physics 
was argued in the same context given above. They introduced the Dirac type Bogoliubov-de 
Gennes Hamiltonian and discussed the group theoretical aspects of the Cooper pair, the 
mean fields. They also discussed the response of superconductors to circularly polarized 
light [22]. P. Strange discussed the gap equation based on the works of Capelle-Gross [26]. 
Therefore, we wish to extend their theory to a more general form, and discuss various aspects 
of the relativistic superconductivity. 

Recently, the relativistic superfluid theory is developed in quantum hadrodynamics ( 
the theory for the system of protons, neutrons and mesons ) [27,28]. Also recently, in the 
quark matter in quantum chromodynamics ( QCD ), the possibility of the formation of 
the Cooper pair ( the color-superconductivity ) is discussed [29,30]. These theories treat 
the superconductivity and/or superfluidity in relativistic many-body systems. The neutron 
stars [31] may become a subject of these theories. The superconductivity in the high-density 
plasma should also be treated in relativistic theory, though we do not find any of this attempt 
in literature. 

With these tendencies of recent theoretical physics, and considerations given above, we 
investigate the BCS and generalized BCS theories in relativistic quantum field theory in this 
paper. Needless to say, the most important approach in the theory of superconductivity is 
the BCS and generalized BCS schemes. The essential idea of the BCS theory is constructed 
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on an effective model Lagrangian, like the Nambu-Jona-Lasinio model or the Gross-Neveu 
model [6]. In our work, we introduce the relativistic BCS Lagrangian, and study the dy- 
namical C/(l)-gauge symmetry breaking. The attractive interaction should arise from some 
physical processes in a matter, as discussed in the theory of the mechanism of the origin of the 
superconductivity. Here we obey the attitude of the generalized BCS theory, and construct 
the BCS and generalized BCS theory in the relativistic framework for studying various sys- 
tems. Since we select the gauge freedom as U(l), our main subject is many-electron-positron 
systems. We try to keep a close relation with condensed matter theory. Namely, we select 
the simplest formalism and intend to study a general theory for the relativistic supercon- 
ductivity. But our formalism can also be applied to the color-superconductivity in QCD. In 
this work, we are mainly interested in the methodological aspect. We want to present how 
the relativistic theory of the superconductivity is constructed and what characteristics it 
has. We intend to construct an unified approach about the superconductivity in the systems 
under strong relativistic effects, from solid state physics to the nuclear and particle physics. 
Due to the standing on the attitude of the BCS and generalized BCS theory, we can realize 
this purpose. 

This paper is organized as follows. In Sec. [TI], the Gor'kov formalism [32] is introduced 
from a Lagragian as a starting point. We regard that this Lagrangian gives nonrelativistic 



BCS theory in its nonrelativistic limit. In Sec. |T|, we perform the group theoretical con- 
sideration of the mean fields in complete manner. In Sec. [IV], especially for treating the 
spin triplet pairing, the generalized BCS theory is constructed on the basis of the previous 
results. In Sec. [V], we solve our Gor'kov equations which have several different types of 



the mean fields, and obtain the Green's functions in analytical forms. In Sec. VI, we give 



the gap equations for various states and discuss their features. Finally in Sec. |V1I| , we 
give summary of this work, and also give a plan about coming works in near future. The 
numerical part will be published as part II of this study. 
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II. RELATIVISTIC GOR'KOV EQUATION 



The purpose of this section is to introduce the basis of the theory of the relativistic 
BCS superconductivity by using the field theoretical Green's function method. For this 
procedure, following the nonrelativistic theory, we use the canonical formalism. In particular, 
we select the Gor'kov formalism [32], because it is the most fundamental formalism in 
superconductivity, and also because we want to obtain the Green's function which increases 
the ability of the applications. Another purpose is to construct the theory which will be 
refered to in the generalized BCS formalism. 

First, we introduce the following Lagrangian as our starting point: 



where the first term is the kinetic term of the electromagnetic field, F^ u is the antisymmetric 
electromagnetic tensor, and is the electromagnetic potential. The second term is the 
kinetic term of the Dirac field, the third is the Dirac mass term, ip and ip are the 4-component 
bispinors describing the Dirac fields. We select 7 " 1 " = 7 , 7^ = — 7 l (i = 1,2,3). We use the 
metric convention as = diag(l, —1, —1, —1) and hence 7^7^ + 7^7^ = 2g fiV . The main 
object of our theory is electron-positron gas. Since generally we treat the system in which 
the numbers of electrons and positrons are different ( finite density ), we have to add the 
next term to our Lagrangian, 



as the fifth term. This term is the background to neutralize and stabilize the system. J M is 
a classical current. This term, however, is not used to describe the physics hereafter. The 
fourth term of the Lagrangian is the four-body contact interaction at the same space-time 
point. It will describe attractive interaction between fermions with the coupling constant 
go > 0. This term can be called the relativistic BCS interaction, and we take here the 



C(x) 



--F^(x)F^(x) + ^(x)i7 M (^ - ieA^(x))ip(x) - mi^(x)ip(x) 
+ 9 ^(x)^x)) 2 , 



(1) 



nx)A,(x), 



(2) 
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simplest form. We regard this Lagrangian involves the usual BCS theory in its nonrelativistic 
limit. To make it clear, we here take the standard representation: 

(A 



(3) 



Here, is the large component while x is the small component. The small component 
becomes negligible at the nonrelativistic limit, and therefore the interaction ^ (ipip) 2 be- 
comes ^(0V) 2 - From the view point of the condensed matter physics, this Lagrangian is 
not completly microscopic, such as the electron-phonon theory, but should be regarded as 
a phenomenological effective Lagrangian. However, if once we introduce this Lagrangian, 
the dynamics of the system is completely determined from it. This theory is, of course, un- 
renormalizable, and we have to introduce a cut-off. This Lagrangian itself has symmetries of 
Poincare invariance, [/(l)-gauge invariance, charge-conjugation invariance, spatial inversion 
and time-reversal invariance. 

Now we obtain the field equations from the Lagrangian by the action principle: 

^ iYid^- ieA^ip - mip + g (ipip)ip, (4) 
n = ^_ 8 dC 

= ~i(d^ + ieA ll )ipY - + g {ip^)ip, (5) 

and the Hamiltonian becomes 

H = J d 3 x$(-i^ ■ V + m - e 7 M M - 7 V)V> - y J d 3 x{^) 2 . (6) 

Here the second term will give the attractive interaction for g > 0. Because we consider 
a many-body problem, we introduce the chemical potential \x in the above Hamiltonian. It 
descibes the finite density at /x ^ as the conjugate of the particle number minus antiparticle 
number. Throughout this study, we treat /i as a parameter introduced from the outside of 
the system. Thus we set /j, — 6f ( Fermi level ). Now, to obtain the Gor'kov equation, we 
use the equations of motion ( Eqs. (4) and (5) ) of the following form: 
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«7 M (^ ~~ ieA ii)^ -m^ + g (ipip)ip + yrfi) = 0, 



0\T JT 



(7) 
(8) 



where T means the transposition, and the second equation is obtained from (5) by taking 
the transposition. 

Now we employ the canonical quantization about the Dirac fields to switch to the quan- 
tized fields. Namely, we demand the anticommutation relations 



{^ a {xQ,x)^\(x ,y)} = 5 a/3 5 {3) (x - y), 
{^(x ,a;),^(xo,i/)} = {ipi(xo,x),^l(xo,y)} = 0. 



(9) 
(10) 



Here a and (5 are the spinor indices. We do not quantize but treat it as a classical field. 

Next, we make some preparations to study the theory of superconductivity. First, we 
introduce various propagators. We use the 8-component Nambu notation [33]: 



^ if)(x) 
ip (x) 



, V(x) = ^(x),ii) T (x)). 



(11) 



Here if) and if) are in the Heisenberg representation. The definition of the generalized one- 
particle propagator is 

G(x,y) = -i{T*(x)*(y)) 



' -i(T$ a (x)$fj(y)) -i(T$ a (x)$(y)) ^ 
-i(Tif) a (x)if)p(y)} -i(Tif) a (x)ipT(y)} J 



1 S F (x,y) a/3 -iF(x,y) a p^ 



V 



(12) 



-iF(x,y) aP -S F (y,x) Pa 

This is a 8x8 matrix. T means the time-ordered product, and (• • •) means the sum over the 
expectation value between states coupled through the superconducting pair-correlation. Sp 
is the Feynman propagator for quasiparticle, — iF and — iF are the anomalous propagators. 
Next, we obtain the equations of motion for the propagators (12), in a similar way as the 
nonrelativistic case. We employ the Gor'kov factorization in (7) and (8), taking account 
of only the superconducting pair-correlation by introducing the mean-field approximation. 



Then we obtain the relativistically generalized Gor'kov equation written down as an 8x8 

matrix equation, Li(x)G Gor ' kov (x, y) = lS^(x — y), namely, 
/ . 



S F {x,y) -iF(x,y) * 



V 



K -iF(x,y) -S F (y,x) T J 



(13) 



i^(dn - ieA^) - m + 7 // A (a;) 

A(x) il^id^ + ieA^) + m - 7 0T /z ) 

' 5^{x-y) ^ 
v ^{x-y) j 
A{x) and A(x) are 4x4 matrix mean fields, so called order parameters. The definitions are 

A(x , x) af3 = g F(x£, x; x , x) a(3 = ^(VU^o , x)^(x , x)), (14) 

A(a; ,a;) a/ 3 = #0^0 > x 'i x o, x) a p = S'o^a^o » aO^/jO&o, x ))- ( 15 ) 

This gives the self-consistency conditions. Eq. (13) is determined self-consistently. Needless 
to say, we can obtain infinite order series of the Dyson-type equation from (13). In the 
superconducting state ( or under the critical temperature ), the mean fields have finite 
values, and we stand on the theory which does not conserve the particle number and charge; 
we attain different representation by non-perturbative way. In general, the mean field clearly 
violates the Lorentz symmetry, as well as the gauge symmetry. In other words, the mean 
field involves quantities other than the scalar. In the case of "superconductivity" theory 
in particle physics, the object is the vacuum: They treat only the scalar value (ipip) 7^ 
( the Hartree fields ) and they discuss the chiral symmetry breaking [2,5]. These points 
are different from our theory. We would like to mention the relation between our Gor'kov 
equation (13) and the relativistic Bogoliubov-de Gennes equation which was given in ref. [18]. 
If we explicitly write the Coulomb potential of nuclei for of (13), our theory becomes 
essentially the same as the relativistic Bogoliubov-de Gennes theory. 

We will also obtain the Fourier transform of the Gor'kov equation for the homogeneous 
system. We set the external field A^ = 0, and then obtain it in the matrix form: 

^ ( S F (k) -iF(k) ^ ^ 1 n\ 



ft - m + 7°/i A 

A f + m- 7 0T /i 



[-iF(k) -S F (-k) T J 



1 



(16) 



Here ft T means the transpose of l/k. The self-consistency condition now becomes 



d 4 p - 



(17) 
(18) 



Because of the homogeneity, the mean field has only the internal degrees of freedom. In the 
nonrelativistic BCS theory, the mean field has no internal degree of freedom. In the case of 
the relativistic theory, there is a possibility to obtain much more complicated states. 

The finite-temperature theory of the Matsubara formalism can be obtained in the same 
way. We introduce imaginary time r = it. The temperature Green's functions are defined 
as 

g(x, y ) = -{TM^(y)) 



' -(T T ^ a (x)$ p (y)) -(T T Mx)^(y)}^ 



\ 



(19) 



-F{x,y)aj3 S(y : x)p a J 
where (• • •) means the statistical average. From the equations of motion of the temperature 
Green's functions, the Gor'kov equation becomes 

A(x) 



^ -7°(gf - A*) +il h dk -m 



S(x,y) af3 -F(x,y) a!3 ^ 



V 



A(x) -7° T (|: + A*) + il kT d k + m 



V 



-F{x,y) a p -S(y,x)p a 



) 



\ 



(20) 



5^(x-y) 

5^{x-y) ) 

Here the definition of the mean fields are the simple extension of those for the zero temper- 
ature: 



A(r, x) afi = g F(r + , x] r, x) = g (^ a (T + , x)$£(t, x)), 



_T 



A(T,x) a p = g J r (r + ,x;T,x) = g (ip a (T + , x)^(r, x)). 



(21) 
(22) 



Fourier transform for the homogeneous system are also obtained, and the Gor'kov equation 
in the matrix form becomes 



^ 7°(iu; n + //) — 7 • k — m 



S(u n ,k) -f(w w fc) 



V 



7° T (ic<; n — ji) — 7 T ■ + m 



\ 



-F(uj n ,k) -S(-u n , -k) T J 



(23) 



Here (3 = l/k B T ( the Boltzmann constant ), w n = (2n + is a fermion discrete 

frequency. Solving the finite-temperature Gor'kov equation, we will obtain the solutions of 
the same form, except that we need to substitute ko — > iu n . The solutions of the Gor'kov 
equation for homogeneous system are discussed in Sec. V. 

III. GROUP THEORETICAL CONSIDERATION OF THE MEAN FIELDS 

We proceed the group theoretical consideration so as to treat the mean fields more easily. 
From the anticommutation relation for Fermi fields, the mean fields obey 



A(x) 



af3 



-A(x)pa, A(x) a p = -A(x)/3 a . 



(24) 



Note that A is not the Hermitian conjugate of A: 



A(x) = 9o(^ T (xMx)) = g ol OT (^(x)^(x)h° = 7°A t (x) 7 °, 



(25) 



but rather, 7°A is the Hermitian conjugate of 7°A. On the other hand, from the gauge 
transformation: 



ip(x) -> tp'(x) = e ia{x) ip(x), ip(x) -> tp'(x) = i){x)e- iai - x \ 



(26) 



the one-body propagator is transformed as 

e iM*)-<*{v))S F (x,y) e i{a{x)+a{y)) {-i)F(x,y) 

e - i{a (x)+a(y)) (-i)p( X , y) e -^)-^v)) (-)S F {y, x) T 



S F (x,y) -iF(x,y) ] 



-iF(x,y) -S F (y,x) J 



J 



\ 



(27) 
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and the mean fields are transformed as 

A(z) Q/3 - e 2 ^A(xU, A(x) a(S - e- 2M(l »A(x)^. (28) 
Under the spatial inversion, 

ip(x ,x) 7V(zo, -a;), i>(x Q ,x) ^ 4>(x , -03)7°, (29) 
the mean fields are transformed as 

A(x , a) = (W T ) Z 7° W T >7° = 7°A(x , -a ; )7 , (30) 
A(x , a) = (^> - 7°$ W = 7°A(x , -x)^ . (31) 

Under the proper Lorentz transformation ( e^ v = —e UfM , <j^ u = 7"] ), 

4>'(x') = Sty(z), ^'(a/) = ^(x)S'- 1 , (32) 
5 = exp(-^0, 5" 1 = 7 ^ t 7°, (33) 

the mean fields are transformed as 

A'(x') = <^V T ') = {S^ifF) = SA{x)S T 
= (l-^)A(x)(l-J v Kf) 

= A(x) - ^„K", A^V^V, (34) 
A' (*') = A(x) + ^e^V^VA^O. (35) 

Here we use C~ lr y^C = — 7^ T , C~ x a^ v C = —cr^ uT and charge conjugation matrix C = 
27 2 7°. Then, taking account of the parity, we find we may take scalar as 7 1 7 3 = — 75C ( 
75 = «7°7 1 7 2 7 3 ) and pseudo scalar as 7°7 2 . We expand the 4x4 matrix mean field into the 
16-dimensional complete set of 7- matrices [19,34]: 

A a/3 = {A s l + A^ + A T ^ + A^ 7 57 M + A p i l5 }(- l5 C). (36) 

In this expansion, we take a notation that S denotes the scalar, V denotes the vector, T 
denotes the 2-rank antisymmetric tensor, A denotes the axial vector and P denotes the 
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pseudo scalar. It is clear from our derivation that each set of this expansion is linearly 
transformed into themselves. This structure of the mean fields in the relativistic theory is 
one of the essential differences from that of the nonrelativistic theory, as discussed in Ref. 
[19]. 

Here we have to mention about the matrix structure of each term given in (36). The 
scalar, pseudo scalar and vector are antisymmetric with respect to transposition, while other 
cases are symmetric. These matrix structures have to be related to the Pauli principle, as 
given in (24). If we want to treat the cases of the axial vector and 2-rank antisymmetric 
tensor, we have to extend our treatment, and the mean fields have to possess finite angular 
momentum with odd-parity. This can be treated by introducing the generalized BCS scheme, 
which will be given in the next section. 

Under the charge conjugation, if) and iJj are transformed as 

V> ^ C^ T , $ ±> -^ T C~\ (37) 

and hence A and A are transformed as 

A = ^ C(i) T i))C- 1 = CAC-\ (38) 

A = $ T $) ^ C^ T )C- X = CAC- 1 . (39) 

Under the time reversal, T = i'j 1 ^ 3 , with the relation 

i>(x ) ^ TiP(-x ), $(x ) ^ ${-x )T, (40) 

together with the rule of taking the complex conjugate about c- numbers, we obtain 

A = (^(xoW T (x )) ^ -T(^(-x )^ T (-x ))*T = — TA*T, (41) 
A = ($ T (x )j>(x )) ^ -T^ T (-x )^(-x )}*T = -T(A)*T. (42) 

This transformation is related to the concept of "unitary" of the mean field [14] . Throughout 
this study, we will treat only unitary pairing state. 

Chiral transformation can also be treated in the same way, 
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i/j(x) -> il)'{x) = e^'^^ix), rl>(x) -> ^'0*0 = V>(zK 7Ba(x) , 



(43) 



then ( = 75 ) we obtain 



1750(0;) 



A' = (ip'ix)^ 1 (x)} = e l ^ x > (ip(x)ilj 1 (x))e l ^ a{x > = e l ^ a(x > A(x)e 
A'(x) = ^ T \x)^' \x)) = e^ a{x) (x)^(x))e iTMx) = e^ a{x) A{x)e^ 5a(x) . 



(44) 
(45) 



Thus we will find that, by using the expanded form of the mean field ( Eq. (36) ), the 
vector and the axial vector are invariant under the chiral transformation. In short, in case 
of the scalar and pseudo scalar, A' = e 2 * 750 ^ A holds, in case of the vector and axial vector, 
A^ = A M holds, and in case of the antisymmetric tensor, A'^ v = e 2ll5a ^ A^ v holds. 

IV. THE GENERALIZED BCS THEORY 



In this section, to treat especially the spin triplet Cooper pair states, we investigate the 
method of the generalized BCS approach in the relativistic theory. Like the nonrelativistic 
theory, our starting point of studying the generalized BCS theory is a Gor'kov equation with 
generalized pairing scheme. We introduce next generalized Gor'kov equations: 

(46) 
(47) 
(48) 
(49) 



(i 7 "^ -m + i°fi) x S F (x, y)-ij d 4 zV(x, z)F(x, z)F(z, y) = 6& (x-y), 
(iy^ -m + -y°fi) x F(x, y)-if d 4 zV(x, z)F(x, z)S F (y, z) = 0, 
(iT" 3 ^ +m- l 0T fi) x F(x, y)+ij d 4 zV(x, z)F(x, z)S F (z, y) T = 0, 
(il^d» + m--f r ii) x (-)S F (y,x) T -i J d 4 zV(x, z)F(x, z)F(z,y) = 5 (4 \x - y). 



In the Fourier transform in a matrix form, they are expressed as 

^ ( S F (k) -iF{k) \ ( 1 o\ 



^-m + 7°/i A(k) 
[ A(k) f + m-^n ) y 
where the definition of the mean fields are 



-iF(k) S F (-k) rj 



J 



v 01 ; 



A(k) = J 
A(k) = 



d 4 p 



(2n) 4 
d 4 p 



V(k,p)F(p), 
4 V(k,p)F(p). 



(50) 



(51) 
(52) 
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Here V(k,p) is an effective attractive interaction which is assumed to give rise to super- 
conducting states. Here we consider V(k,p) as a scalar function. We can assume it as a 
Feynman propagator for a scalar field. The definite form of the function V(k,p) is related 
to the symmetry of the Cooper pair; different symmetries correspond to the different forms 
of V{k,p). Since the matrix structure of (50) is same as (16), the Green's function obtained 
by solving (50) gives the same matrix structure. It is also clear that this can be extended 
to the finite-temperature Matsubara formalism. 

Next, we decompose the interaction and the mean fields into each channel. When we 
denote the total momentum of a Cooper pair as P M = (Po, P), P 2 = Pq—P 2 , and when P 2 > 
0, we can always stand on a rest frame which satisfies the condition P = 0. We treat our 
problem under the condition P > 0, P = 0. Our mean fields have the translation invariance. 
At that time the Wigner little group of the Poincare group becomes 0(3) rotation [34]. It 
should be enough to use the irreducible representation of 0(3), i.e. the three dimensional 
spherical harmonics, for the channel decomposition. It is also related to the fact that our 
theory will mainly treat the finite-temperature Matsubara formalism, which selects a specific 
time-coordinate. In this formalism, the theory allows only 0(3) rotational symmetry. Under 
such consideration, we assume that the interaction depends only on the angle between k 
and p, expanding it by using the addition theorem, 

v(k, P ) = ^^i + iMkoAk^poApm^Y^ip) 

l mi 

->EE 4 *fc(P). (53) 
i m 

where we also introduce the weak coupling approximation Vi(ko, \k\;po, \p\) ~ gi [35]. This 
approximation neglects the retardation of the interaction and also neglects the dependence 
on the magnitude of each momentum. We may call this as "anisotropic contact" interaction. 
If one of these channels is attractive, the Fermi sea can be unstable and it can become a 
superconducting state. 

Later we will find that the gap equations which can have nontrivial solutions are those 
for the scalar A 5 , the Oth component of vector Ajf, the space-like components of axial 
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vector A A , the axial- vector-like components of 2-rank antisymmetric tensor A^. We have 
to discuss here only these mean fields. To satisfy the Pauli principle, for the case of spin 
singlet pairing, we select the parity-even part of the expansion of the interaction (53), 

V^(k,p) = g + g 2 f^ Y 2m (k)Y* m (P) + • • • ■ (54) 

m=2 

Here the first term corresponds to the interaction g (^ilj) 2 /2. The mean field is also expanded 
in the same way, 

-2 

A 1 = A + £ A 2m r 2m (fc) + • • • . (55) 

On the other hand, for the spin triplet pairing, we have to select the parity-odd part of the 
expansion. As for the interaction, we have 

V^(k,p) =5iE Y lm (k)Y* m (p) + 93 E Y Zm {k)Y; m {p) + ■■■. (56) 

m=l m=3 

We need some consideration for the treatment of the spin triplet mean fields. Because of 
the relativistic nature of the theory, spin and orbital degrees of freedom can not be treated 
separately. The spin projection to the space-fixed z-axis is not conserved. Now the good 
quantum numbers are the total angular momentum j, its projection rrij and helicity A. 
We will treat the spin triplet Cooper pair state as an axial vector boson, with even internal 
parity. The basis for the expansion of the triplet mean field then becomes the helicity states. 
We define the space-fixed coodinate as the k x k y k z -system, while the body-fixed coordinate 
moving with the boson as the {r/C-system. We choose the direction k along the £-axis. The 
helicity is the projection to this axis, A = —1, +1 (transverse), A = (longitudinal). Here we 
take the phase convention after the textbook of Landau and Lifshitz [36]. Then the basis 
becomes 

Wfc) = i^^^eWpgiCfc), (57) 
where e^ A ^ is the spherical unit vector of the ^r/C-system, 

e( ±1 ) = =F^(|±^), e<°>=iC, (58) 
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As for the Wigner rotation matrix, we determine the relation between the Eulerian angles 
and spherical angles of the k x k y k z -system as 

0,7) = Z>2i(«M,0) = e^d^(9). (59) 
The orthonormal relation is given from the equations, 

e^)* . e (A 2 ) = <W (60) 

J\in9d9 j* ' d4>V%»(k)V^(k) = (61) 

as 

/ sin 9d9 d(f)ip* imiXl (k) ■ ij j2m2 x 2 (k) (62) 

" J 

About the inversion a = (ft ^ <fi + n , (3 = 9 ^ n — 9, 7 ^ 7r — 7, each function is transformed 

as 

= (-iy- x v«L(-k), (63) 

eW(fc) = (-^W^H-fc). (64) 
Thus the parity eigenstates become 



^mo(k) = ii-y^±l e (°)2>S£(fc), (65) 
V$A|(*) = V-'^^i^V^k) + e ( - 1) P°L(fc)), (66) 
d| W = i j - 1 ] l^(e^(k) - e^V^ m (k)). (67) 

Each parity is given as; ip jm0 : : ("i)^ 1 , ^jm|A| : ( _1 ) j - Thus the odd-parity 

states becomes ^000 f° r the monopole (j = 0), V'imi f° r the dipole (j = 1), ^ 2m o and V4mi f° r 
the quadrepole, etc. Since in this work, we treat only the unitary states of the mean fields, 
and (A-cr)(A* -a) = A ■ A* +i(A x A*) -a, we have the unitary condition A x A* = [14]. 
The "longitudinal" ipj m o is itself unitary, but as for the "transverse" V'jmi, we have to add 
the complex conjugate to itself. This means we have to take a linear combination with the 
time-reversal state in order to obtain a unitary state. Using the relations 
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= (-i) J - m ^L, 
C } r = (-i)^ m+ Vj:L, (68) 



we obtain the following bases that we use for the expansion: 



i^jmo, (69) 



^[^i + (-ir m ^Li], (70) 
1 

V2 l 



-^[^-(-ir^tli]- (7i) 



In the "longitudinal" basis, the direction of the mean field vector coincides with k, while 
in the "transversal" basis, it is orthogonal to k. Finally the spin triplet mean fields are 
expanded in the following form: 



A(fc) = 5>^;m,o(fc), (72) 

= £ A« _L(^w l( fe) + (-iy-^^ mjl (k)), (73) 

A(fc) = E A&-^&> x(fe) - (-ly-^'^iW), (74) 

where we have to select only the parity-odd part of the expansion, for satisfing the Pauli 
principle. Here we also take the so-called weak coupling approximation Aj m (k , \k\) m 
A jmj [35]. 

V. THE SOLUTIONS FOR THE HOMOGENEOUS SYSTEM 

Now using the definition C _1 7 M C = — 7 mT and the notation k = (k + fc), = 
(A; — /i,k), our Gor'kov equation becomes 

(j& - m)S F (k) - iAF(k) = 1, (75) 

— m)F(k) — iASF{—k) T = 0, (76) 

— C _1 (^ — m)CF(k) + iASp{k) = 0, (77) 

-C -1 ^ - m)C(-)5 F (-A;) T - iAF(k) = I. (78) 
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First, S F (k) and F(k) satisfy the relation 



Then we obtain 



F(k) = ^ l i-^-CAS F (k). (79) 
k — m z 



(jk - m)S F (k) + AC' 1 1 + m CAS F {k) = 1. (80) 

k 2 — m 2 

T 



On the other hand, S F (—k) T and F(k) satisfy 



F(k) = i M^- 2 AS F (-kf, (81) 



then 



(ft - m){CS F {-k) T C- 1 ) + CA f + m AC-^CSFi-kfC- 1 ) = 1. (82) 

A; 2 — m 2 

Therefore, our problem turns out to solve the two 4x4 matrix equations (80) and (82). Un- 
fortunately, it is difficult to solve these equations completely because of the matrix structure 
of the Dirac operator, particularly in the analytical form. We have to solve the equations 
assuming the type of the mean field that might be realized. We use the expanded form of 
the mean field, 

A = {A s + A p i lb + A V + A^ 757 ° 

+A y • 7 + A A ■ 757 + Af y) • * 7 °7 + Af A) • 75 7°7}(-7 5 C), (83) 

and put each type of the mean field separately into the Gor'kov equations. In this way, we 
can treat each cases of the mean fields in the same manner. In particular, the case of the 
vector mean fields, under the unitary assumption, can be treated completely in the same 
way as the scalar mean field case. We write down the expanded form of the mean field ( 
Eq. (83) ) as, 



AC- 1 = -Q2A a .r a }j B , (84) 

a=l 

CA = - 75 {E A : • r °}- (85) 



a=l 
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Going back to (80) and (82), puting one of 8 types of mean fields A a , under the unitary 
assumption, (A a -T a )(A* -T") = (A* -r a )(A a -r a ) = ±(A a -A*) (where plus sign corresponds 
to the scalar, 0th component of vector, space-like components of axial vector and axial- 
vector-like components of 2-rank antisymmetric tensor, while minus sign corresponds to the 
other cases ), and also using ^ 5 T a = ±r a 7 5 ( plus sign corresponds to the scalar, pseudo 
scalar and antisymmetric tensor, while minus sign corresponds to the vector and axial vector 
), we obtain 

(jf - m)S F (k) - (A a ■ nl^(A: ■ nS F (k) = 1, (86) 

k z — m z 

$ - m^Spi-kfC- 1 - (A: • r a )J^- 2 (A a ■ nCS^-kfC- 1 = 1. (87) 



k 2 — m 



Solving this equations, we obtain 



S F (k) = j^{(k 2 - m 2 )§ + m) - (A a ■ ml + m)(A* a • F 1 )}, 

-iF(k) = (±)^y7 5 {(^ - m)(A a • r a )(# + m) - (±)(A a • A* a ){A a ■ T a )}C, 

-iF(k) = ^C-HW ~ m)(A* a • r°)(jf + m) - (±)(A a ■ A:)(A: ■ P 1 )}, 

-S F (-k) T = -^f-'iCk 2 - m 2 )(lt + m) - (A* a ■ r»)§ + m)(A a • r°)}C, 
D(k) = (k 2 - m 2 )(k 2 - m 2 ) 

- {j{A a ■ v a )UK ■ n + (A a • nhK ■ ni - 2m 2 (±)(A a . a:)} 

+ (Aa • A:) 2 . (88) 

As for the sign for F(k), plus corresponds to the scalar, pseudo scalar and antisymmetric ten- 
sor, while minus corresponds to the vector and axial vector. All the eight Green's functions 
we have obtained reflect the features of each type of the mean fields, both in the numerator 
and denominator. We can confirm that this solution satisfies the Gor'kov equation. We can 
see that when we take the limit A a — > 0, the solution gives the Green's functions of normal 
state. 

We have solved the Gor'kov equations under the assumption of specific types of the mean 
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fields in order to make them tractable. We have to check the existence of coupling among 
different types of the mean fields. For this purpose, we take the gap equations 

A(*) = / j^V(k,p)F(p). (90) 

We substitute the expanded form of A ( Eq. (83) ) into the left hand side, substitute the 
solutions for F(p) obtained previously into the right hand side, and take the trace of both 
side. We find couplings between A s and Aq, A p and A„ , A v and A™, A a and AT A y We 
should say that the Green's functions obtained above are approximate ones by neglecting 
these couplings, and also the gap equations obtained by using these Green's functions contain 
these approximations. 

In the end of this section, we investigate the factorization of the denominators of the 
Green's functions. We discuss some features of quasiparticle dispersion. In all of the cases, 
the denominator given in (88) becomes the following form: 

D(k) = (k 2 - m 2 )(k 2 - m 2 ) - 2(A a • A* a )(k ■ k + d) + (A a ■ A* a ) 2 . (91) 

This is 2nd order in k 2 , thus easily factorize it as 

D{k) = (k - E + )(k + E + )(k - E„)(k + EJ). (92) 

D(k) is fourth order in k , but the Gor'kov equation itself is 8x8 matrix. Like the non- 
relativistic 1 S'-BCS theory, the quasiparticle spectra are doubly degenerate. We will find 
that the condition for doubly degenerated quasiparticle spectra is the unitary of the mean 
fields, which implies the time reversal symmetry of the mean fields [14]. E + is the branch 
for the quasiparticle coming from the positive energy solution, while E_ is the branch for 
quasiparticle coming from the negative energy solution. The expression for the pole becomes 

E 2 ± = k 2 + m 2 + fi 2 + (A a ■ A:) =f yV(fc 2 + m 2) + 2 (A a • A* a )(m^ + d), (93) 

where d = —m 2 for the scalar, d = +m 2 for the pseudo scalar, d = 2k 2 — m 2 for the 0th 
component of vector, d = 2\k- A y | 2 /(A y - A v *)+m 2 for the space-like components of vector, 
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d = 2k 2 + m 2 for the Oth component of axial vector, d = 2\k ■ A A \ 2 / (A A ■ A A *) — m 2 for the 
space-like components of axial vector, d = —2\k ■ AT v ^\ 2 /(AT v ^ • A£M + 2k 2 + m 2 for the 
vector-like components of 2-rank antisymmetric tensor and d — —2\k- Aj A ^\ 2 / (Aj A ^ ■ Aj A ^) + 
2k 2 — m 2 for the axial-vector-like components of the 2-rank antisymmetric tensor. About 
the d, the sign before m 2 is coming from the parity of the mean fields which is determined by 
7-matrices, minus for even parity and plus for odd parity. It holds that in the limit m — > 0, 
the quasiparticle spectra for the scalar and pseudo scalar, Oth component of the vector and 
axial vector, space-like components of the vector and axial vector, vector-like components 
and axial- vector-like components of the 2-rank antisymmetric tensor, coinside. We write the 
quasiparticle dispersion explicitly. For the scalar obtaied in Ref. [18], 

E 2 ± = (^k 2 + m 2 Tfi) 2 +\A s \ 2 . (94) 

In the pseudo scalar case, 

E 2 ± = k 2 + m 2 + fi 2 + \A a \ 2 t 2\/(k 2 + m 2 )fi 2 + m 2 \A p \ 2 , (95) 

and these spectra coinside at m — ► 0, El = (\k\ =f /i) 2 + |A| 2 . In the case of the Oth 
component of the vector pairing, 

El = k 2 + m 2 + /i 2 + |A^| 2 t 2\f(k 2 + m 2 )fi 2 + \A^\ 2 k 2 , (96) 

and in the limit m — > 0, El = (\k\^f^J jj? + (Aq'I 2 ) 2 . In the case of the space-like components 
of the vector pairing, we obtain 

El = k 2 + m 2 + /i 2 + (A v ■ A v *) t 2\/(k 2 + m 2 )^ 2 + m 2 (A^ • A y *) + jfe • A^| 2 . 

(97) 

In the case of the Oth component of the axial vector pairing, 

El = (\/k 2 + m 2 t VW|A^| 2 ) 2 , (98) 

and in the limit m — > 0, the cases of the Oth component of the vector and axial vector 
coincide with each other. The case for the space-like components of the axial vector pairing, 
we obtain 
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E 2 ± = k 2 + m 2 + fi 2 + (A A ■ A A *) t 2^J(k 2 + m 2 )fi 2 + \k ■ A A \ 2 , 

(99) 

and for m — > 0, this also coinsides with the case of the space-like components of the vector 
pairing. About the vector-like components of the 2-rank antisymmetric tensor, we obtain 

El = k 2 + m 2 + /i 2 + ( Af y) • Af; } ) T ^(k 2 + m 2 )^ 2 + Aj v) • A^) - \k • Af v) P, 

(100) 

and the axial-vector-like components of the 2-rank antisymmetric tensor becomes 

El = k 2 + m 2 + /i 2 + (Af A) • Af; } ) =F 2 v /(fc 2 + m y + k 2 (Af A) • A^) - |fc • Af A) p. 

(101) 

These spectra also coinside with each other in the limit m — > 0. 

The dispersion of the quasiparticle includes a gap topology and it reflects orbital part 
of the mean field. It is clear that the complexity of the dispersion in relativistic case is the 
result of the matrix structure of the Dirac operator. We can make an assertion that these 
situations manifest the fact that in the relativistic theory, we cannot decouple the spin and 
orbital degrees of freedom. The dispersion and also the gap topology reflect the density of 
states and then affect the thermodynamics like the specific heat. They influence also the 
collective modes. Spin part of the mean field gives spin symmetry of the Cooper pair, and 
it is reflected to the spin susceptibility or the Knight shift. In all the cases given above, we 
obtain El = {\Jk 2 + m 2 =F /i) 2 in the limit A a -> 0. 

VI. THE GAP EQUATIONS 

In this section, we will obtain gap equations by using our Green's functions obtained 
previously. We use the finite-temperature Matsubara formalism. The gap function deter- 
mines the thermodynamic properties of the superconductivity as direct consequences. The 
gap equation itself is the stationary condition of the free energy [10]. 
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Under the Pauli principle, and with the condition gi > 0, we will find that we have 
nontrivial solutions for the mean fields in the case of the scalar and 0th component of vector 
for even-parity interaction V^ e \ the space-like components of axial vector, and axial- vector- 
like components of 2-rank antisymmetric tensor for odd-parity interaction V^°\ Other cases 
give the gap equation of the form "1 = negative value", thus under gi > 0, there is no 
nontrivial solution. Thus in the interaction g (ipi/j) 2 /2 with g > 0, it can give nontrivial 
solutions only in the case of the scalar pairing and 0th component of vector pairing. 

First, we deel with the interaction go^ijj) 2 /2. The self-consistency condition becomes 

A = *? £ *>• (102) 

1 r rl 3 k 

£ = S oX^/^K,*). (103) 

Substitute T for scalar or 0th component of vector into (102), deal with the 7-matrices and 
take trace of both sides. Hereafter we only deal with the nontrivial case A a 7^ 0. In the 
denominator D(u n , k) = {ioj n — E + ){iuj n + E + )(iu n — E^){iu n + E_) ( substitute k — > iu n 
in (92) for the Matsubara formalism ), we perform the partial fraction decomposition and 
discrete frequency summation. For the scalar case, we obtain the gap equation of the same 
form as derived by Strange [26]: 

l = ^/ A 7 ^(J-tanh^ + + J-tanh^A (104) 
2 J- a (2tt) 3 V2£ + 2 + 2£_ 2 /' K J 

E± = \j (\/ k 2 + m 2 t l^) 2 + |A 5 | 2 . (105) 

Here we introduce the cut-off for the momentum integration. Of course, this cut-off scheme 
can not be a covariant form. In the high-energy physics, there are problems that a different 
cut-off scheme gives different physical quantities [4,5]. But we do not treat these problems 
here. Eqs. (104) and (105) give simple extension of the nonrelativistic 1 S'-BCS theory. The 
second term of the integrand in (104) is the contribution of the quasiparticle coming from 
the negative energy solution. 

In the case of the 0th component of the vector pairing, the gap equation is given in the 
following form: 
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2 M 2 ") H ^(fc 2 + mV + | A V| 2fc 2 ; 2£+ 2 
+{l + ^= } — — tanh-ffV (106) 



£ ± = ]/ k 2 + m 2 + fi 2 + \A%\ 2 t 2^(k 2 + m 2 )fi 2 + \A%\ 2 k 2 . (107) 

On the other hand, in the case of the pseudo scalar pairing and 0th component of the 
axial vector pairing, the gap equations become 

g f d 3 k /( m 2 -> 1 ft 



- 1 



( jl }— -tanh-ff 4 

VI //, 2 , ...ON. .9 , I A PI9„_9i 2K 9 ^ 



2 7 (2tt) 3 VI v /( fc 2 + m 2 )/U 2 + | A P|2 m2 J2^ 

+ ; m " = |— tanh^EA (108) 

1 ^/(fc 2 + mV + | A P|2 m 2/2E_ 2 ^ 



and 



£ ± = \k 2 + m 2 + y 2 + |A P | 2 t 2 V(fe + ™V + |A p | 2 m 2 , (109) 



l = ^/-^( / 1 tanh^ + + . 1 tanh^A (110) 

2i(27T)3V 2 ^2 + | A A|2 2 2^2+|A^ 2 ^' 

Vfe 2 + m 2 =F^ 2 + |A^| 2 . (Ill) 



We find that the right hand side of each equations are positive. Hence, we do not obtain 
nontrivial solutions in these cases. Due to the same reason, the cases for the space-like 
components of the vector pairing and the vector-like components of the 2-rank antisymmetric 
tensor pairing also have no nontrivial solutions for gi > 0. 

Next, we treat the spin triplet states by using the generalized BCS formalism. Adopting 
the Green's functions of the space-like components of the axial vector or the axial-vector- 
like components of the 2-rank antisymmetric tensor for the self-consistency condition, and 
manipulating the 7-matrices, we obtain a 3-vector gap equation: 



x{A(p)(p 2 Tp 2 -m 2 - fI 2 - A(p) • A») ± 2p(A(p) • p)}. (112) 
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Here the upper sign corresponds to the case of the axial vector while the lower sign corre- 
sponds to the case of the axial-vector-like components of the 2-rank antisymmetric tensor. 
We take the expansions of the interaction and mean fields as discussed previously. Through- 
out our study, we take only I = 1 p-wave interaction into account. For the mean fields, we 
treat j = 0, 1, 2 states. We can consider coupled states with different j, rrij, A, but hereafter 
we only treat definite values of j, rrij, A, because of simplicity. In the general form, 

d 3 p _ 1 „ ,~ _ 1 



x{Aj^ jmjA (p)(p^p 2 - m 2 - /i 2 - |V; mjA (p) • ^ jmj x(P)) 
±p(A^ (^a(p)-p))}, (113) 

we multiply both sides with ip*j m .\(k) from the left, and perform angular integration with 
respect to the spherical angle of k. After a lengthy manipulation, we obtain 18 gap equations, 
and all of them have the forms of the following 2-type gap equations, either Type-(i) or 
Type-(ii): 
Type-(i) 

1 = / i^Mi 1 - I " 2 w ) dr teIlh I E + 

" • (p 2 • • &|A$ V(^) 



+ ( X + / P m ) 4" tanh f ^ }' (1M) 

V ^(p 2 WV 2 + &|A$/p 2 c(0) J 

£± = + ™ 2 + /i 2 + &|A$/c(0) =F 2^(p* + m 2 )/i2 + & |A^.| Vc(^). (115) 
Type-(ii) 

l = a0i / t^c(^) f — tanh -£+ + — tanh ^£_Y (116) 

y i (2tt) 3 v ; V2£ + 2 + 2£_ 2 /' v ; 

£ ± = yJ(y/fP + m*Tri 2 + b\AW ,1*0(0). (117) 

In the above equations, a, b and c(0) depend on specific ipj mj x- All the 18 gap equations are 
listed in table I. In this table, we also list the node structure of the gap functions. 
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About -^ooo cases, there is no angular dependence both in the eigenvalues and gap equa- 
tions. The condensation is for j = I + s = pairs, ^ooo state of corresponds to the BW 
( Balian-Werthamer ) state in our theory [9,35]. The gap is isotropic like the nonrelativistic 
BW state. The structure of this gap equation is the same as the A s given above ( Eqs. 
(104) and (105) ). Vioi state of A A can be regard as the ABM ( Anderson-Brinkman-Morel 
) state in our theory [8,35]. In this case, the spin vector is orthogonal to the vector of orbital 
motion. The angular dependence of this gap in the gap equation is the same as Y\±\, like 
the nonrelativistic ABM. It has nodes at two points 9 = 0, ir. There is no state of j = 1 
which has a line node structure. It is interesting that, in spite of the spin triplet pairing, 
the structures of the gap nodes of j = 2 states are like that of d-wave, Y 2m . This feature 
arises from our relativistic treatment, i.e. the use of the helicity state. 

In the spin-triplet pairing states discussed above, the quasiparticle dispersion and the gap 
equation take simple forms like the case of A s except the angular dependences of the gaps 
( Type-(ii) ), for the axial vector when the orientation of the mean fields is perpendicular 
to k ( helicity ±1, the "transversal" ), and for the 2-rank antisymmetric tensor when the 
direction of the mean fields is parallel with k ( helicity 0, the "longitudinal"). Other cases, 
the eigenvalue and the gap equation have complicated structure in Type-(i), and they are 
similar to that of the Ajf. These are clear from the form of the functions, (99), (101) and 
(112). Thus in the relativistic theory, the difference in the helicity ( it gives the orientation 
of the mean field vector ) gives rise to a quite different form for both the dispersion and gap 
equation. It might be interesting if the relative orientation between the mean field and k 
gives large effects to the solutions of the gap equations, thermal properties or response of 
the system to the external fields. These gap equations will be studied numerically in part 
II of this paper. 
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VII. SUMMARY 



In this paper, we have performed the investigation of BCS and generalized BCS super- 
conductivity in relativistic quantum field theory. We have introduced the Gor'kov equation, 
and given the group theoretical considerations about the superconducting mean fields in 
complete manner. Especially to treat the spin triplet Cooper pairs, we have investigated the 
generalized BCS formalism in our theory. We have solved the Gor'kov equations completely 
under the assumption to have the specific types of the mean fields with unitary condition. 
We have constructed the gap equations in the various states and discussed their details. 
Throughout this paper, we have discussed various characteristic features of the theory. 

As a next work, we will present the results of solving the gap equations numerically. 
We will also give some results of the thermodynamics, or response to the external fields. 
Preparation for the presentation of these results is now in progress, and will be published 
as part II of this study. 
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TABLES 



TABLE I. The coefficients in the gap equations of the A A and AjT^ pairings, (i) and (ii) 
correspond to the use of the gap equations of Type-(i) and Type-(ii), in (114,115) and (116,117), 
respectively, a, b and c{9) in this table are the coefficients in these equations. 







type 










state 




AT 
A (A) 


0, 





c{V) 


node 


^000 


l 1 ) 


( n ) 


l 

2 


J_ 

47T 


I 


no node 


^200 


W 


(n) 


1 

4 


5 
16tt 


(3 COS P — 1) 


a _l l 

COS0 = ± 73 


vbn n 


d) 
\ x j 


(ii) 

v. ./ 


3 
2 


15 

87T 


sin 2 9 cos 2 


= o, | )7r 


^220 


(i) 


(ii) 


3 
8 


15 
32tt 


sin 4 # 


6" = 0, 7T 




(ii) 


(i) 


3 
4 


3 

8tt 


sin 2 fl 


= 0,7T 




(ii) 


(i) 


3 
8 


3 
16tt 


cos 2 9 + 1 


no node 




(ii) 


(i) 


9 
4 


15 
8tt 


sin 2 6* cos 2 9 


= 0,§,7T 




(") 


(i) 


3 
8 


5 
16tt 


4cos 4 #- 3cos 2 fl + 1 


no node 




(ii) 


(i) 


3 
8 


5 
16tt 


1 - cos 4 


6> = 0,7T 
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